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Abstract. We investigate the optical conductivity in the Mott insulating phase of the one-dimensional 
extended Hubbard model with alternating hopping terms (dimerization) at quarter band filling. Optical 
spectra are calculated for the various parameter regimes using the dynamical density-matrix renormal- 
ization group method. The study of limiting cases allows us to explain the various structures found nu- 
merically in the optical conductivity of this model. Our calculations show that the dimerization and the 
nearest- neighbor repulsion determine the main features of the spectrum. The on-site repulsion plays only 
a secondary role. We discuss the consequences of our results for the theory of the optical conductivity in 
the Bechgaard salts. 
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1 Introduction 

The eleetronic properties of quasi-one-dimensional charge- 
transfer salts have been intensively investigated in recent 
years ^ElEl- An important example of such compounds 
is the family of Bechgaard salts {TM)2X, where TM is 
the organic molecule TMTSF (tetramethyltetraselenaful- 
valene) or TMTTF (tetramethyltetrathiafulvalene), and 
X denotes an anion such as ClO^, PFg", Br~, etc. These 
materials have highly anisotropic structures and proper- 
ties. Therefore, it is believed that above an energy scale 
of a few meV their electronic properties can be described 
in first approximation by one-dimensional models. 

The one-dimensional Hubbard model with alternating 
hopping integrals (dimerization) and a quarter-filled band 
has been proposed to describe various properties of the 
Bechgaard salts [311|6l l x||8] | 9] . in particular their unusual 
optical spectrum |l()lllllll2lll8lll4j . In this model the for- 
mation of a Mott insulating ground state is due to the in- 
terplay of the Coulomb interaction between electrons and 
the lattice dimerization. However, the relevance of this 
approach for Bechgaard salts has remained controversial. 
Although we know the generic features of the low-energy 
optical spectrum in one-dimensional Mott insulators [151 
I16II17| . the optical conductivity of the quarter-filled Hub- 
bard model with dimerization has not been determined 
accurately yet. Thus, no direct comparison with the ex- 
perimental spectrum observed in Bechgaard salts has been 
possible. 

The recent development of the dynamical density ma- 
trix renormalization-group (DDMRG) method |15II18| al- 



lows us to calculate the dynamical properties of low di- 
mensional correlated electron systems with an accuracy 
comparable to exact diagonalizations but for much larger 
system sizes. Here, we apply the DDMRG method to the 
calculation of the optical conductivity in the quarter-filled 
dimerized Hubbard model. The model and method are 
briefly introduced in the next section. Then, we present 
our results in Sec. 01 Finally, we discuss the consequences 
of our results for the theory of the Bechgaard salts in 
SecH 



2 Model and Method 

The one-dimensional extended and dimerized Hubbard 
model is defined by the Hamiltonian 



H 



E (4 

odd l;(7 
*2 E {^.A+l.a 



(1) 



even L-.a 



It describes fermions with spin cr =1 , j which can hop be- 
tween neighboring sites representing the highest occupied 
molecular orbital (HOMO) of each TM molecule. There 
are three electrons in the HOMOs of each pair (TM)2, so 
that the band made of the HOMOs is three-quarter filled 
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in terms of electrons or quarter filled in terms of holes. We 
use the hole representation and keep the number of par- 
ticles N such that we have a density p = N/L = 1/2 for 
an even number of lattice sites L. The operator c^^ (q.o-) 
creates (annihilates) a hole with spin a at site I. The hole 
density operator is fii^a = c^^ci^a and fii ~ fii,^ + rijj is 
the total number of holes at site I. The hopping integrals 
ti > t2 > give rise to a single-particle dispersion 



e(fc) = ±JA^ sin^(fc) + cos^{k) 



(2) 



with a total band width At = 2ti+2t2 and a (dimerization) 
gap 2Z\ = 2{ti — ^2)- The Coulomb repulsion is mimicked 
by a local Hubbard interaction U, and a nearest-neighbor 
interaction V. The physically relevant parameter regime 
for Bechgaard salts is U > 2V > 0. In Tabled we show 
some values of the model parameters ti,t2,U, and V which 
have been proposed [S||n|[7||19) to describe various (TM)2X 
salts. 

We use open boundary conditions since density ma- 
trix renormalization group (DMRG) algorithms are most 
efficient for this type of boundary |2()II21) . In open chains 
it is important to use the correct form of the (non-local) 
Coulomb interaction between electrons in the Hamilto- 
nian (QJ. Neglecting the average density (i.e., setting p — 
in Eq. ^ results in complicated edge effects in the excita- 
tion spectrum such as the existence of low-energy excita- 
tions localized at the chain ends. 

Two mechanisms can induce an insulating ground state 
in this model at quarter filling jHUHl- First, the Umklapp 
scattering in an effectively half-filled band [—2t < e{k) < 
—A] due to the dimerization can lead to a Mott insulat- 
ing state accompanied by a ikp bond order wave (BOW), 
where hp = t^p/'^ = 7r/4 is the Fermi vector. Second, for 
large enough parameters U and V the Umklapp scattering 
in the quarter-filled band [—2t < e(fc) < 2t, neglecting the 
dimerization gap 2A] can drive the ground state into an 
insulating phase with a spontaneously broken symmetry: 
a Akp charge density wave (CDW) P22n23j . In the fam- 
ily of Bechgaard salts, TMTSF compounds are believed 
to be realizations of one-dimensional Mott insulators ^Hj 
while TMTTF compounds are considered to be charge 
ordered [5| like in a CDW state. For realistic parameters 
(see Table^, however, the system described by the Hamil- 
tonian is a Mott insulator [H]. Therefore, we will in- 
vestigate the optical conductivity in the Mott insulating 
phase only. Since we use open boundary conditions, we ob- 
serve 2fci?-B0W and 2kF- and 4fci^'-CDW fluctuations in- 
duced by the chain ends (Friedel charge oscillations) in the 
ground state. For all the parameters U,V,ti,t2 discussed 
in this work, however, the ground state has no long-range 
order or broken symmetry but the 4fci?-B0W induced by 
the alternating hopping terms ti ^t2- 

To determine the ground state properties and to ob- 
tain some information about excited states of the Hamilto- 
nian |^ we use a standard DMRG technique |20II2H . For 
instance, the Mott gap (also called single-particle charge 
gap or charge transfer gap) 



Table 1. Model parameters (in meV) for various Bechgaard 
salts from Refs. |SI6l'rli^| . 









ti 


U 


V 


(TMTSF)2PF6 


(Ref. El) 


250 


225 


1250 





(TMTSF)2C104 


(Ref. |71) 


290 


260 


1450 


210 


(TMTTF)2PF6 


(Refs. laini) 


135 


95 


945 


380 



can be calculated from the ground state energies Eo{Nh) 
for Nfi holes in the system ^Cj. 

The linear optical absorption is proportional to the real 
part (Ti (w) of the optical conductivity, which is related 
to the imaginary part of the current-current correlation 
function by 



CTi(w > 0) = -ilm {tPolJ ^ 



J\^o). (4) 



Here, \ipo) is the ground state of the Hamiltonian H, Eq 
is the ground state energy, and 77 — > 0+. Assuming that 
the sites are equidistant, the current operator J is 



odd l:a 



-it2 {4,a^l+l,a - cj+i. 



(5) 



even l\a 



With these definitions the optical conductivity ax{uj) is 
given in units of e'^a/h, where 2a is the lattice constant 
and e the charge of a hole. The frequency uj is given in 
units of t/h. 

In an open chain the optical conductivity is also re- 
lated to the imaginary part of the dipole-dipole correlation 
function 

Mlo) = ^lui {^o\D- ^ ^^D\^o) , (6) 

^ Eq + UJ + n] — H 



where the dipole operator is 



D 



1=1 



1=1 



L + l 



ni. 



(7) 



Ec = EoiN + 1) + Eo{N - 1) - 2Eo{N) 



(3) 



One can apply the DDMRG method to chains of 
finite size L to compute the optical conductivity ai(uj) 
with a finite broadening r] > 0. Thus, DDMRG yields 
the convolution of cri(w) with a Lorentzian of width rj or 
the quantity defined by Eq. ^ or Eq. [S] for a finite ij (see 
Ref. for more details and the advantages of the vari- 
ous approaches). The properties of the optical spectrum in 
the thermodynamic limit L ^ 00 can be determined us- 
ing a finite-size-scaling analysis flHj with an appropriate 
broadening 

V{L) ^ II L. (8) 

This approach has already been successfully used to study 
the optical properties of simple one-dimensional Mott in- 
sulators (i.e, in the extended Hubbard model at half fill- 
ing) |15ll24ll25j . In particular, a quantitative description 
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has been achieved |26) for the experimental low-energy op- 
tical conductivity spectrum in the quasi-one-dimensional 
compound SrCu02. 

Often one can use deconvolution techniques to com- 
pute a smooth spectrum without broadening from the 
numerical DDMRG data for finite rj and finite system 
size |27II28| . In this work we use a standard linear reg- 
ularization method for the inverse problem I 29J to decon- 
volve DDMRG spectra. The deconvolution usually yields 
a very accurate description of the spectrum in the ther- 
modynamic limit if it does not possess any sharp feature 
(i.e., on a scale smaller than the broadening ry used in the 
DDMRG calculation). Therefore, the broadening r] used in 
the DDMRG calculation sets the resolution of a spectrum 
obtained through a deconvolution. 

An extension of the DDMRG method ^Hl can be used 
to compute the excited states of the Hamiltonian which 
contribute to the optical spectrum Q). We have used this 
method to determine the optical gap (i.e., the excitation 
energy uji of the lowest eigenstate with a finite matrix 
element {ipi\J\ipo)) in finite chains more accurately. 

Up to m = 320 density-matrix eigcnstates have been 
kept per block in DDMRG calculations and up to m = 768 
in ground state DMRG calculations. Truncation errors are 
negligible for all results presented here. Thus, the accu- 
racy of our calculations is mostly limited by the finite 
broadening or resolution 77 ~ 1/L imposed by finite sys- 
tem lengths. 



3 Results 

In the absence of the non-local electron-electron interac- 
tion {V = 0), the properties of the model (Q) depend on 
the parameters U and A only. (The average hopping term 
t just fixes the energy scale.) First, we will investigate 
three limiting cases ^ for which the main features in the 
optical conductivity cri(ti;) can be easily understood: (i) 
the large-dimerization limit t2 <C ti(=> A 2t), U < 4ti, 
(n) the strong-couphng hmit U :s> ti > t2, and (in) the 
weak-coupling limit U ^2 < ti. Then we will discuss 
how the optical spectrum changes when the parameters U 
and A are varied between the limiting cases. Finally, we 
will consider the effects of the nearest-neighbor repulsion 
V. 



3.1 Large Dimerization 

In the dimer limit A ~ 2t (t2/ii 0) the system is 
made of independent dimers (i.e., pairs of nearest-neighbor 
sites) 0. The system eigenstates are products of the dimer 
eigenstates (i.e., the eigenstates of a two-site Hubbard 
model). In the ground state at quarter filling each dimer 
is occupied by exactly one hole, which is localized on 
that dimer. The current operator ((S)) does not couple the 
dimers. Therefore, only intra-dimer excitations can con- 
tribute to the optical conductivity. (They correspond to 
transitions from the bonding orbital to the anti-bonding 




2 4 6 8 10 



as 

Fig. 1. Optical conductivity ai{ui) in the large-dimerization 
limit (Zi = 1.64t) for a strong effective coupling {U — 3.64f « 
2042) with a broadening rj — 0.05t {L = 128 sites). Inset: same 
data on a logarithmic scale. 



orbital of the dimer in the limit U = 0.) It can be shown 
that <7i{uj) consists of a single Dirac (5-pcak at w = 2ti ~ 
At = 2A for any U > 0. Note, that moving one hole from 
a dimer to another one (inter-dimer excitations) yields 
eigenstates of the system with an excitation energy which 
can be lower than 2A and thus the Mott gap Ec is lower 
than the optical gap wi = 2A in that special case. (For 
instance, Ec vanishes as U goes to zero but wi = 2 A > 0.) 

We now discuss the optical excitations for a finite inter- 
dimer hopping t2 ^ ti and U < iti. (For larger U/ti the 
spectrum is better understood starting from the strong- 
coupling limit U ^ ti, which is discussed in the next sec- 
tion.) For small but finite t2 the dimer eigenstates become 
hybridized and build bands of delocalizcd electronic states 
with a bandwidth oc t2- Thus, in the optical spectrum the 
(5-peak at w = 2ti is replaced by an narrow absorption 
band ('intra-dimer' band) of width cx ^2 around lo = 2ti 
[approximately for 2A = 2{ti ~t2)<uj<4:t = 2(ti -1-^2)]- 
Figures and [3 show the optical conductivity calculated 
with DDMRG for A = IMt (=» ^2 = O.lSt and ^1/^2 ~ 
10) and two different couplings U = 3.64t and U = 0.546i, 
respectively. The 'intra-dimer' band contains a substan- 
tial part of the optical weight and is clearly visible as the 
strong feature at3.4<w<4.2in Figs. ^ and |2| 

The current operator now couples nearest neighbor 
dimers with a term cx t2. Thus, for finite t2 inter-dimer 
excitations also contribute to the optical spectrum. At 
high energy (oj > 2A) these excitations give rise to two 
small peaks around oj = 2ti and oj = U + 2ti. These 
features are much weaker than the 'intra-dimer' band be- 
tween 2A and At as their optical weight is of the order of 
t\/ti and t\/{U + 2ti), respectively. Nevertheless, the first 
peak is clearly visible on the top of the 'intra-dimer' band 
in Fig. n and the second peak in the inset of that figure 
at w ~ 7.5. We note that these excitations correspond to 
moving one particle from the bonding orbital of a dimer to 
the anti-bonding orbital of another dimer. In particular, 
the inter-dimer excitation which appears in the middle of 
the 'intra-dimer' band at cj = 2ti corresponds to the for- 
mation of a triplet state on the second dimer. Clearly this 
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Fig. 2. Optical conductivity ai{u!) in the large-dimerization 
limit {A — 1.64f) for a weak effective coupling {U — 0.546t « 
3t2) with a broadening rj = 0.2t {L = 64 sites). Inset: high- 
resolution and expanded view of cti {to) in the low-energy region 
UJ < 0.2t. DDMRG results (circles) for r? = 0.0128f (L = 200 
sites) and field-theoretical result I15| (line) for a gap Ec = 
0.049t and the same broadening r;. 



optical excitation involves both charge and spin degrees 
of freedom. 

The low-energy spectrum (uj < 2A) is more interest- 
ing. In the large-dimerization limit the model ^ can be 
mapped onto a half-filled Hubbard chain with effective 
parameters teS = ^2/2 and Ucs = U/2 for U small com- 
pared to 4ti Consequently, the low-energy spectrum is 
given by the optical conductivity ai{iu) of the half-filled 
Hubbard model, which is known |15| . For instance, for 
A = 1.64t and U = 3.64i the effective interaction is strong 
C^cff/^cff = U/t2 « 20. Accordingly, the shape of the low- 
energy band (w < 2) in Fig. ^ is similar to the semi- 
elliptic absorption band centered around uj = Uos = 1.82t 
found in the strong-coupling limit of the half-filled Hub- 
bard model [T31ISn] . Also, the optical weight in this struc- 
ture is of the order of tlg/Uce = t\/U and thus much 
lower than in the 'intra-dimer' band. For A = 1.64t and 
U = 0.54t, however, the effective interaction is relatively 
weak ?7off/tcff = U /t2 ~ 3.3 which corresponds to a small 
Mott gap Ec ~ 0.049t ~ 0.54tcff- In that case, the optical 
weight in the low-energy region lo < 2A is significantly 
larger than for a strong effective coupling and comparable 
to the weight at higher energy, as seen in Fig. [3 Moreover, 
the low-energy spectrum calculated with DDMRG (shown 
with a higher resolution in the inset of Fig.|21) agrees very 
well with the field-theoretical prediction in the limit of a 
small Mott gap [T5] . 



3.2 Strong Coupling 

We now turn to the strong-coupling limit U 3> > ^2- 
In this regime the dimerization opens gaps of 2A in the 
lower and upper Hubbard bands ^ with band widths 
At. At quarter-filling low-energy elementary charge excita- 
tions (holons in the lower Hubbard band) have the same 
dispersion ^ as electrons in a half-filled band (Pcicrls) 




Fig. 3. Reduced optical conductivity ujai{uL>) in the strong- 
coupling limit {U = 40t) calculated with A = 0.6t and a broad- 
ening 77 = 0.2t (L = 64 sites). The dashed line is the Peierls 
insulator spectrum for the same values of A and rj. Inset: ex- 
panded view of the high-energy spectrum. 



insulator. In particular, there is a Mott gap Ec ~ 2A. Ne- 
glecting the contribution of the spin degrees of freedom 
to the matrix elements {ipn\J\'>Po)j where {ipn) is an ex- 
cited state, one expects ^H] that the optical conductivity 
(Ji{uj < U) is similar to that of a band (Pcicrls) insula- 
tor EH 



(71 (w) 



(2Z\)2(4t)2 



4tjV[^^ - (2Zi)2][(4t)2 - w2 



(9) 



In Fig.Olwe compare this analytical result with our numer- 
ical DDMRG data for U = 40i and A = 0.6t. Both spectra 
have been broadened with a Lorentzian of width r] = 0.2t 
to facilitate the comparison. The agreement is excellent 
but for a small shift ~ /U, which can be attributed to 
the finite value of U used in the numerical calculations. 
At high energy u! > U there is also a weak absorption 
band with a total spectral weight oc t'^ /U corresponding 
to charge excitations from the lower to the upper Hubbard 
band (see the inset of Fig. |2l . 

The predictions of the strong-coupling theory remain 
qualitatively valid for relatively weak couplings U. For in- 
stance, the main features of the spectrum jnj are clearly 
visible in the optical conductivity shown in Fig. 01 which 
has been calculated with DDMRG for U = 5ti and ti/t2 = 
2 (corresponding to U/t = 20/3 and A/t = 2/3). There is 
also a weak absorption band not described by Eq. © at 
high energy to > U (see inset of Fig. 0J as in the strong 
coupling limit. Note, however, that there are clear quan- 
titative differences. For instance, for these parameters we 
have found a Mott gap Ec = 0.53t (in full agreement with 
Ref. ^), which is clearly smaller than the strong-coupling 
result 2 A « l.'iit. 

A close inspection of the DDMRG spectra reveals a 
weak structure at low frequency 2A < w < 2t which is 
not explained by the simple theory predicting the spectral 
form ij^. This feature is barely visible around w = 3 on 
the scale of Fig. 13 but can be clearly seen as a small bump 
around lu = 2.8 in Fig. ^ To understand this deviation 
from Eq. it is helpful to analyze the spectrum as a 
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Fig. 4. Optical conductivity ai{uj) for A/t = 2/3, U/t = 20/3, 
and a broadening ry = O.lt {L = 128 sites). Inset: expanded 
view of the high-energy spectrum. 

function of the dimcrization parameter A for very strong 
coupling U. For A <^ t (i.e., ti w i2) most of the optical 
weight is concentrated in the low-energy singularity at a; = 
Ec = 2 A, which in the limit Z\ — * becomes the Drude 
peak of the metallic ground state. For Z\ « 2i (i.e., ^2/^1 ^ 
1) the optical weight becomes equally distributed between 
both divergences at to = Ec = 2A and co = At. Most of the 
optical weight is concentrated in this narrow band which 
is the counterpart of the 'intra-dimer' band found in the 
large-dimcrization limit (see previous section). Thus, the 
optical spectrum is dominated by a similar structure in 
both the strong-coupling regime [U ^ ti) with A k, 2t 
and the large-dimcrization limit [A » 2t but U < Ati). 
The weak spectral features, however, are quite different in 
both regimes. In particular, there is no optical absorption 
at low energy u; < 2A in the strong-coupling limit (C/ ^ 
^i). The crossover from one regime to the other one is 
quite complicated and will not be discussed here because 
it is not relevant for the (TM)2X salts. 

Nevertheless, comparing the results for large dimcr- 
ization and those for strong coupling we observe that the 
unexplained weak feature in ai(cj) for U ^ ti corresponds 
to the excitation involving both charge and spin degrees 
of freedom at an energy uj = 2ti in the dimer limit. There- 
fore, we conclude that the spin degrees of freedom are re- 
sponsible for the (small) deviation from the simple Peierls 
spectrum © in the large-C/ limit. A similarly small con- 
tribution of the spin degrees of freedom to the optical 
spectrum has already been observed in analytical calcu- 
lations [30] and DDMRG simulations for the half-filled 
Hubbard model [T3] . 

3.3 Weak Coupling 

In the weak-coupling limit t/ <C ^2 < ^i, the low-energy 
sector of the Hamiltonian can be mapped onto a half- 
filled chain with an effective (bare) band width 2t2 and 
an effective long-range electron-electron interaction oc U 
which induces a small Mott gap Ec <^ t2 The low- 
energy properties Ec < co 2t2 of this system should be 
well described by field-theoretical approaches. In partic- 
ular, one expects the optical conductivity to be given by 




Fig. 5. Optical conductivity for A = 0.105t, U = 5.263t, and 
ri = O.lt (L = 128 sites). Inset: same data on a double loga- 
rithmic scale. 



the field-theoretical result for one-dimensional Mott insu- 
lators in the small gap regime jl5lll6j . For weak coupling 
the optical weight must be concentrated at low energy 
uj ^ Ec ^ t, A as one expects that most of the spec- 
tral weight lies in the Drude peak when the system be- 
comes metallic {Ec 0). Therefore, the field theoretical 
approach describes the essential part of the optical spec- 
trum. 

We have performed several DDMRG calculations of 
(Ti(w) in this weak-coupling regime. For instance. Fig. 
shows the optical conductivity calculated for U = 5.263i 
and A ~ 0.105t (corresponding to ^2/^1 = 0.9). Clearly, 
the optical weight is concentrated in a sharp peak at low 
frequency as expected (the long tail at high frequency is 
mostly due to the broadening rj = O.lt). For these param- 
eters the Mott gap is Ec ~ 0.03t in the thermodynamic 
limit and the optical gap converges to the same value as 
seen in Fig. |S1 The broadening of the DDMRG spectrum 
also results in an apparent shift of the peak position [i.e., 
the maximum Wmax of o'i(ti')] to higher frequencies. In the 
limit of an infinite chain and with the scaling (jH)) one finds 
(see Fig. O that tj,„ax approaches a value (0.04t) only 
slightly larger than the Mott gap Ec in agreement with 
the field theory prediction |15lllf>j . As most of the spectral 
weight is concentrated on a scale to Ec comparable or 
smaller than our typical resolution 77, it is not possible to 
make a quantitative comparison between field theory and 
numerical results for the spectral lineshape. Nevertheless, 
our DDMRG results are always qualitatively compatible 
with ficld-thcorctical predictions for the behavior of cti (w) 
at frequencies of the order of the charge gap Ec. 

A field-theoretical analysis |E1 of the optical conduc- 
tivity in one-dimensional Mott insulators predicts that the 
leading asymptotic behavior at high frequency lu ^ Ec is 

aiiio) ^Lo-'', (10) 

where the exponent a > 1 depends on the interaction 
strength. This power-law behavior is a good approxima- 
tion at extremely high frequencies u; ^ lO^Ec only [TH] 
and the field theory approach is only valid for energies 
much smaller than the effective band width ~ 2t. Thus, 
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0.02 0.04 0.06 



1/L 

Fig. 6. Mott gap Ec (circle), optical gap uji (square), and 
position Umax of the conductivity maximum (triangle) as a 
function of the inverse system size for A = 0.105t, U — 5.263t, 
and V = 0. Lines are quadratic fits. 



in the lattice model ^ one can observe such an asymp- 
totic behavior in the limit of small Mott gaps only (i.e., 
for Ec <^ uj ^ t). Moreover, the high-frequency behavior 
can be modified by various processes which are neglected 
in the field-theoretical approach but generate additional 
optical transitions for cu > Ec, such as inter-band tran- 
sitions at uj > 2A or transitions between the lower and 
upper Hubbard band around uj — U . For instance, wc 
have seen in Sec. lH.ll that in the (effective) weak-coupling 
regime of the dimerized limit the low-frequency spectrum 
UJ ^ Ec <^ t can be described with field theory but the 
high-frequency spectrum is dominated by the inter-dimer 
excitations around oj w 2 A (see Fig. Obviously, a 
power-law behavior cannot be observed for Ec uj ^ t 
in that case, fn practice, only the weak-coupling regime of 
the Hamiltonian ^ seems to fulfill both conditions nec- 
essary for the occurrence of the asymptotic power-law be- 
havior of ai{uj): (i) a small Mott gap Ec <^t and (ii) no 
other optical excitation in the relevant frequency range. 

We have found that DDMRG spectra for finite broad- 
ening rj and system size often decay as a power-law at high 
frequency. In the inset of Fig. |S1 one clearly sees such a 
behavior with an exponent « —1.2 for 0.2i < uj < lOt cor- 
responding to TEc < UJ < 330Ec. (In Sec. 18.51 wc will see 
that the high-frequency spectrum is actually dominated 
by another feature explained by a strong-coupling anal- 
ysis.) However, the exponent and the range over which 
the power-law behavior can be observed depend on the 
method used to broaden the spectrum in the DDMRG 
calculation (see Sec. ISJ, the broadening 77, and even the 
system size. Therefore, this power-law decay is probably 
an artifact of our numerical approach. This effect can eas- 
ily be understood if one assumes that most of the optical 
weight is concentrated in a sharp structure at w ~ i?c. The 
broadening of this structure creates a broad tail which de- 
creases asymptotically as r]AiUj~^ , where the exponent (3 
lies between 1 and 3 (depending on the precise broaden- 
ing technique used in the DDMRG simulation) and the 
coefficient Ai is proportional to the total optical weight. 
Obviously, this artificial broad tail hides an asymptotic 




Fig. 7. Optical conductivity ai{ui) for U — 6t and 77 — 0.2t 
{L — 32) for various dimerizations A. 

behavior ai{uj) ~ " for /3 < a. It can also hide 

the asymptotic behavior of <Ji{uj) up to relatively large 
frequencies for (3 > a if the high-frequency spectrum con- 
tains only a small fraction of the total optical weight (i.e., 
A2 -C 11 Ai). We think that this effect is responsible for the 
power-law observed in our DDMRG spectra with a finite 
broadening 77. 

To determine the true asymptotic behavior of (7i{uj) 
we have tried to deconvolve the DDMRG spectra for fi- 
nite 77 in order to obtain spectra for r; = in the thermo- 
dynamic limit |27| . The resulting spectra do not show a 
power-law behavior in any significant range of frequencies. 
Unfortunately, the accuracy of the deconvolved spectra is 
very poor at high frequencies because our deconvolution 
technique (linear regularization approach for an inverse 
problem |29| ) does not work well when the spectrum is 
dominated by sharp structures as in Fig. [S] 

In summary, we have not been able to determine the 
asymptotic behavior of cti (uj) in the weak-coupling regime. 
While some of the raw DDMRG data clearly exhibit a 
power-law behavior at high frequency, we think that this 
is an artifact of our method. We can not confirm (or refute) 
the validity of the field-theory prediction 1)1 for the lat- 
tice model investigated here. Nevertheless, our investi- 
gation leads us to conclude that an asymptotic power-law 
behavior (|10|l can occur only in the weak-coupling regime. 
Moreover, the optical weight at high-frequency (i.e., in the 
asymptotic tail) can only be a small fraction of the total 
optical weight, which is concentrated just above the gap 
Ec. 



3.4 From Small to Large Dimerization 

Although the nature of the optical excitations greatly dif- 
fers in the strong and weak coupling limits, we have found 
that the evolution of the optical spectrum with A is qual- 
itatively similar for all values of J7 > 0. In the strong- 
coupling limit the spectrum is given by Eq. and thus 
the evolution of cri(w) with A can be described accurately. 
When the dimerization is weak {A <^ t), most of the opti- 
cal weight is concentrated just above the first singularity 
at the spectrum onset uj = Ec ^ 2A and the second singu- 
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Fig. 8. Reduced optical conductivity ioai{Lj) calculated with 
rj = 0.2t (L — 64) for A = 1.64t and various values of U. Inset: 
expanded view of the high-frequency conductivity spectrum. 

larity a,t u; = At carries very little weight. As A increases, 
the optical weight is progressively transferee! from the low- 
energy singularity to the high-energy one until the spectral 
weight becomes equally distributed between both singu- 
larities as one reaches the large dimerization limit A — > 2t. 
Simultaneously, the first singularity moves to higher en- 
ergy as A increases and ultimately merges with the second 
(fixed) one as A reaches 2t. Therefore, we observe both 
a transfer of optical weight from a low-energy structure 
around lo = Ec to sl high-energy structure around lo ~ At 
and a shift of the low-energy structure toward higher ex- 
citation energies as A increases. 

Away from the strong-coupling limit Eq. ^ is not an 
accurate description of the optical spectrum. Nevertheless, 
our DDMRG calculations show a qualitatively similar evo- 
lution of the spectrum as a function of A for all values of 
U that we have analyzed (i.e., down to U = t). For in- 
stance, Fig. 13 shows the optical conductivity calculated 
with DDMRG for U = Qt and several values of Z\. We 
clearly see both the optical weight transfer from the low- 
energy peak to the high-energy structure and the shift of 
the low-energy peak toward higher energy as A increases. 
We note, however, that the low-energy peak is close to the 
Mott gap Ec for small A only. For larger A the position 
of the first peak moves away from Ec (at least when U is 
not too large) contrary to the strong-coupling result ©. 
As a result, the dominant features in tri [uj] can lie well 
above the Mott gap Ec as already shown for the large- 
dimerization limit in Sec. 13.11 The high-energy structure 
always lies at an energy close to At for all U and A but 
its weight can become so small that it is no longer visible 
such as in the weak-coupling limit (see Fig. |3J). 

3.5 From Weak to Strong Coupling 

For a given dimerization A < 2t the strength U of the 
local Coulomb interaction significantly modifies the Mott 
gap Ec, which is equal to the optical gap for A < 2t and 
V = 0. However, it has little effect on the main features of 
the optical spectrum and just modifies the fine structure. 
For instance, Fig.|Hlshows the reduced optical conductivity 
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Fig. 9. Reduced optical conductivity ujai{Lo) calculated with 
rj — 0.2t (L — 64) for A — O.lObt and various values of U. Inset: 
expanded view of the high-frequency conductivity spectrum as 
a function of lj — (7. 

calculated with DDMRG for A = l.GAt and various values 
of U. For U < At the optical spectrum is well explained 
by the analysis of the large-dimerization limit fSec. I^.ljl . 
It consists of a strong structure at a; w 2ti = 3.6At (the 
intra-dimer band) and weaker features at uj fa Ec and 
w ~ 2ti + U (see Fig. Pl. When U increases, the gap 
becomes larger and correspondingly we observe a progres- 
sive shift of the low-frequency weak structure toward the 
strong peak in Fig. |H1 Simultaneously, the high-frequency 
weak feature moves to higher energies in agreement with 
the relation uj « 2ti + U (see the inset of Fig.lSjl- However, 
the strong dominant structure remains largely unaffected 
by the variation of the coupling U. 

Nevertheless, the analysis of the optical conductivity 
for changing U yields an interesting result. As discussed 
in Sec. 13.31 in the weak-coupling approach the local in- 
teraction term is responsible for an (effective) interaction 
oc U which splits the (effectively) half-filled band [defined 
by the lower part of the single-particle dispersion (0)] in 
two (effective) Hubbard bands separated by a gap Ec- Ex- 
citations from the lower to the upper effective Hubbard 
bands significantly contribute to the optical spectrum on 
the energy scale uj = Ec- In the strong-coupling approach 
the local interaction term splits the full band defined by 
the dispersion Q in two (full) Hubbard bands separated 
by a gap ~ U. In that case, excitations from the lower 
to the upper full Hubbard bands contribute to the opti- 
cal spectrum around uj = U. Our calculations show that 
both features can be seen (for instance, in Figs. ^ and |HJl 
for a given value of U. Therefore, we conclude that, at 
least for some parameters {A,U) of the model QJ, the 
low-frequency part of the optical spectrum is explained 
by weak-coupling (i.e., field-theoretical) approaches while 
the high-frequency part is explained by a strong-coupling 
analysis. 

As a second example of the optical spectrum evolution 
with U, we show the reduced optical conductivity calcu- 
lated with DDMRG for a small dimerization A = O.lOSt 
and various interaction strengths U in Fig. |51 For U = t 
the system is in the weak-coupling limit (discussed in 
Sec. l3.3f) and most of the spectral weight is concentrated in 
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Fig. 10. Optical conductivity cri(tj) for A = 0.1054, U = 
5.263t, rj = Q.2t {L = 64), and various nearest-neighbor inter- 
actions V. 

a peak close to ui ^ Ec- A second weaker feature is visible 
at about w = 4i (see the inset of Fig. (Sj) and corresponds 
to transitions from the bottom to the top of the single- 
particle band For larger U the optical weight remains 
concentrated in the low-energy peak at uj ~ Ec- This peak 
moves to slightly higher energies and appears to broaden 
because the energy scale set by the gap Ec increases with 
U until it reaches 2 A for [/ ^ oo as discussed in Sec. 13.21 
but its shape is not significantly changed by the variation 
of U. The feature at w ~ 4i disappears for U > 4t but 
another weak feature becomes visible around lo U for 
strong enough coupling U (see the inset of Fig.EJ. Again 
this corresponds to optical excitations from the lower to 
the upper (full) Hubbard bands. 

We note that for U = 5.263t this contribution to the 
optical spectrum is already clearly visible around uj — U = 
At in the inset of Fig. This result contradicts the appar- 
ent asymptotic power-law decrease (IIUII discussed previ- 
ously for the same parameters (see Fig. EJ. This discrep- 
ancy is due to the different broadening methods and data 
representation used for Figs. [S] and |5| It confirms that the 
asymptotic power law found in some of our spectra for 
weak couplings are probably an artifact of the broaden- 
ing used in the DDMRG calculation. This also illustrates 
how difficult it is to observe the field-theoretical predic- 
tions (|10|l for the asymptotic behavior of ai{uj) in the lat- 
tice model because optical transitions neglected in the 
field theory approach significantly contributes to the high- 
frequency spectrum. 

In summary, our calculations for the model |^ with 
V = show that the distribution of the optical weight 
is essentially determined by the dimerization amplitude 
A. For a fixed A only the fine structure of the optical 
spectrum and the energy scale set by the gap Ec depend 
significantly on U . 



3.6 Nearest- Neighbor Interaction 

Neglecting the long-range part of the Coulomb interaction 
is difficult to justify in an insulator. In this section we 
consider the effects of a nearest-neighbor repulsion V > 
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Fig. 11. Mott gap i?c (circle), optical gap lo\ (square), and 
position Wmax of the spectrum maximum (cross) as a function 
of the inverse system size for (a) A — 0.105t, U = 5.263t, and 

V = 2.105t, and (b) A = 0.353t, U = 8.2S5t, V = 3.294i. Lines 
are quadratic fits. 

in the Hamiltonian (This term mimics the long-range 
part of the Coulomb interaction.) For large enough V the 
nature of the ground state of changes from a Mott 
insulator to a CDW insulator [Sj. Here we discuss only 
the optical conductivity in the Mott insulating phase. 

A previous DMRG investigation [7| of the model Q 
has shown that the charge gap Ec increases with V in 
the Mott insulating phase. Our calculations confirm this 
result. For V = (and A < 2t) the optical gap uji is 
equal to the Mott gap Ec in the thermodynamic limit 
(see Fig. for an example) . This gap marks the onset of 
an excitation continuum of unbound pairs of elementary 
charged excitations (for instance, holon-antiholon pairs in 
the weak-coupling picture |16| ) which are responsible for 
the lowest absorption band in the optical conductivity. 
We have found that the optical gap remains equal to the 
Mott gap for non-zero but small nearest-neighbor repul- 
sion V. Thus, the low-energy spectrum still corresponds 
to unbound pairs of charged excitations. 

Figure ^| illustrates two main effects of the nearest- 
neighbor repulsion on the optical spectrum for small V. 
First, the low- frequency peak above Ec is shifted to higher 
frequency as V increases because Ec also increases. Sec- 
ondly, the total spectral weight decreases with increasing 

V as described in Ref. [S]- 

For stronger coupling V, however, we have found that 
the optical gap extrapolates to a smaller value than the 
Mott gap in the limit of an infinite chain. This can be 
seen in Fig. In the first example (for A = O.lOSt, U = 
5.263t, and V = 2.105^) the difference Et ^ lji - Ec 
is small (about O.OAt) but in the second example (for 
A = 0.353t,U = 8.235i, and V = 3.294i) it is large 
{Eb « Q.9t). Moreover, the maximum of the optical spec- 
trum is located at a frequency Wmax which approaches the 
same value as wi in the thermodynamic limit. This is the 
signature of a excitonic peak ((5-peak) at w = wi in the op- 
tical spectrum of the Mott insulator. (See Refs. [T^24II25| 
for a description of the analysis that allows us to identify 
an excitonic peak in a DDMRG spectrum.) The presence 
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of an cxcitonic peak signals a fundamental change in the 
nature of the lowest optical excitation. It is now a neutral 
bound pair (for instance, a bound holon-antiholon pair 
in field theory) called a Mott-Hubbard cxciton j23j. The 
properties of Mott-Hubbard excitons have been investi- 
gated using DDMRG and analytical methods in a previ- 
ous work |24| . The exciton energy oji is smaller than the 
excitation energy of the lowest unbound pair of charged 
excitations in the continuum (which is Ec). The difference 
Eb = oji — Ec is the binding energy of the exciton. There- 
fore, the first example in Fig. ^] corresponds to a weakly 
bound (large) Mott-Hubbard exciton while the second ex- 
ample corresponds to a tightly bound (small) one. 

The exciton generates a (5-peak at to ^ cui < Ec in 
the optical conductivity spectrum. As V increases the ex- 
citon becomes more tightly bound (smaller) and the op- 
tical weight is progressively transfered from the contin- 
uum above Ec to the excitonic peak. In Fig. ^] one sees 
that the optical weight first moves to higher frequency as 
one increases the coupling V form to 1.636i because the 
Mott gap Ec (i.e., the continuum onset) increases. If the 
coupling V is increased further to 3.294t, one finds that 
the spectral weight shifts to lower frequency because it is 
transfered to the excitonic peak at toi ~ 0.5t lying below 
the continuum onset at Ec = lAt (which still increases 
with V). (Note that the gap between the exciton peak 
and the continuum is not visible in Fig.Elfor V = 3.294t 
because of the broadening of the spectrum.) The (weak) 
structure visible at w w 4i in the spectra calculated for 
V ~ rapidly looses weight as V increases. In summary, 
we have found that the nearest-neighbor repulsion V has 
a significant impact on the shape of the optical spectrum 
contrary to the on-site repulsion, but only when it is large 
enough to generate an exciton. 



4 Discussion 

In this section we discuss the implications of our results 
for the theory of the Bechgaard salts. First, we examine 
which values of the model parameters could be appropri- 
ate for (TM)2X salts. Realistic estimates for the hopping 



integrals ti and t2 were proposed more than ten years ago 
on the basis of experimental results and quantum- 
chemistry calculations [l^ ■ Using these estimates Mila |S] 
determined the model parameters U and V from the re- 
duction of the infrared oscillator strength observed ex- 
perimentally in the Bechgaard salts. He found that rela- 
tively large nearest-neighbor repulsions V were necessary 
to explain the reduction of the electron kinetic energy due 
to correlation effects. According to Mila's analysis appro- 
priate model parameters are ^2/^1 = 0.9, U = 5ti, and 
V = 2ti, which correspond to Z\ = O.lOSt, U = 5.263i, 
and V = 2.105t, for (TMTSF)2C104 and iaAi = 0.7, 
U = 7ti, and V = 2.8ti, which correspond to Z\ = 0.353t, 
U = 8.235t, and V = 3.294t, for (TMTTF)2PF6. The op- 
tical spectrum obtained with DDMRG for these param- 
eters are shown in Figs. ^| and El respectively. As ex- 
plained in Sec. 13.61 we have found that the lowest optical 
excitation is an exciton with an energy uji smaller than 
the Mott gap Ec for these parameters (see Fig. 

More recently, Nishimoto et al. [7] fitted the Mott gap 
of the model |^ to the experimental optical gap to deter- 
mine the model parameters. Using the same ratios ^2/^1 
and U /ti as Mila they found that the nearest-neighbor 
repulsions V necessary to reproduce the optical gap data 
were significantly smaller than the values given by the 
reduction of the oscillator strength. According to their 
analysis V = 0.764t for (TMTSF)2C104 and V = 1.636i 
for (TMTTF)2PF6. The optical spectrum obtained with 
DDMRG for these parameters are also shown in Figs. E| 
and^l respectively. In that case we have found that there 
is no exciton and an absorption continuum due to un- 
bound pairs of charged excitations starts &t oj ~ uji ~ Ec- 

The discrepancy between these studies can be under- 
stood. In the (TMTSF)2C104 case the kinetic energy is 
a rather flat function of V for the relevant parameters 
ii, ^2, while the Mott gap increases rapidly with V 0. 
Thus, the uncertainty on Mila's value for V is quite large 
and the value V = 0.764t reported by Nishimoto et al. is 
also compatible with the experimental reduction of the os- 
cillator strength. Therefore, we conclude that for the salt 
(TMTSF)2C104 the nearest-neighbor interaction should 
be close to (though somewhat larger than) the value given 
by Nishimoto et al. and excitons do not play any role in the 
optical excitation spectrum. The appropriate model pa- 
rameters for (TMTSF)2C104 are summarized in Tabled 

In the (TMTTF)2PF6 case, however, the kinetic en- 
ergy is a rather steep function of V for the relevant pa- 
rameters ti,t2,U IHj and the value V = 1.636< found by 
Nishimoto et al. is only compatible with the oscillator 
strength reduction for unrealistically large U. As Nishi- 
moto et al. assumed that the experimental optical gap 
corresponded to the theoretical Mott gap Ec, they ef- 
fectively neglected excitonic contributions to the optical 
spectrum. In particular, their analysis does not take into 
account that the theoretical optical gap cui is significantly 
smaller than the Mott gap Ec for large V when excitons 
occur, as seen in Fig. [rTT b). As a result their analysis 
underestimates the value of V. Therefore, we conclude 
that for (TMTTF)2PF6 the nearest-neighbor interaction 
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should be close to (though somewhat smaller than) the 
value V = 3.294t found by Mila and excitons dominate 
the optical excitation spectrum [at least in the framework 
of the model 1^]. The appropriate model parameters for 
(TMTTF)2PF6 are summarized in Tabled 

We now examine how the main features of the optical 
spectrum in (TM)2X salts can be explained by the dimer- 
ized extended Hubbard model with the parameters deter- 
mined above. Parallel to the stacks of organic molecules 
the optical conductivity of (TMTSF)2X salts has two dis- 
tinct components: a narrow peak at zero frequency (Drude 
peak) with a very small fraction (about 1 %) of the spec- 
tral weight and an absorption band with most of the spec- 
tral weight at finite energy |lllll2lll3lll4j . This second fea- 
ture lies in the mid infrared range above the crossover en- 
ergy above which excitations are effectively confined to a 
single stack and thus can be described by one-dimensional 
models. (Obviously, the zero-energy feature always lies be- 
low such a crossover energy and can only be described in 
the framework of a three-dimensional model.) The finite- 
energy feature is usually interpreted in terms of a Mott 
insulator. When rescaled by the intensity and frequency of 
the spectrum maximum, the optical conductivity of var- 
ious (TMTSF)2X salts exhibit a remarkably similar be- 
havior |13| . In particular, a power law in the frequency 
dependence 

-^,=c(^y'' (11) 

is observed over a decade in frequency 2ci;poak < uj < 

As discussed in Sec. l3.3l our numerical approach is not 
sufficiently accurate to confirm the existence of such an 
asymptotic power-law behavior in the optical spectrum of 
the model iQJ. Nevertheless, our analysis clearly indicates 
that if there is a power-law behavior of (7i{uj) for some 
frequency range in the model the optical weight as- 
sociated with that feature must be extremely small com- 
pared to the optical weight of the peak feature. In the 
experimental spectrum, however, there is substantial op- 
tical weight in the region where the power-law behavior 
is visible. Therefore, we conclude that the universal fea- 
ture jni of the optical spectrum in (TMTSF)2X salts 
cannot be explained within the model 

The well-defined mid infrared structure observed in the 
optical spectrum of the Bechgaard salt (TMTSF)2PF6 is 
difficult to understand in view of the fact that its DC 
conductivity remains metallic down to very low tempera- 
ture. It seems that optical excitations are visible only for 
energies much larger than the energy scale above which 
the system can be seen as metallic (i.e., the Mott gap 
for charge excitations). Therefore, Favand and Mila |S] 
have proposed that the optical gap toi observed in the ab- 
sorption spectrum is much larger than the Mott gap Ec 
because of optical selection rules. Using exact diagonaliza- 
tions of small systems they have argued that such an effect 
occurs in the quarter-filled dimerized Hubbard model |^ 
without the nearest-neighbor repulsion V. Our analysis of 
this model shows that the Mott gap is smaller than the op- 



tical gap only in the dimer case A = 2t (see the discussion 
in Sec. I3.1|l . For all other parameters {ti,t2,U), however, 
we have found that cui = Ec (wi < Ec is also possible for 
V > as shown in Sec. l3.6|l . It can happen that the optical 
weight at the Mott gap is very small while a very strong 
structure is visible in the spectrum at a higher energy. 
For instance, this occurs in the large-dimerization limit as 
seen in Fig. d In a real material with such an absorption 
spectrum the weak low-energy band could easily be over- 
looked, leading to an apparent "optical gap" larger than 
the real gap for charge excitations Ec- For realistic pa- 
rameters, however, we always find that the optical weight 
is very large close to the Mott gap. In particular, for the 
parameters t2 — 0.9ti and U ~ 5ti (A = 0.105^ and 
U = 5.263i) used by Favand and Mila for (TMTSF)2PF6 
(see Tablenj we have found that the Mott gap, the optical 
gap, and the maximum of the spectrum converge to very 
close values in the thermodynamic limit (see Fig. inj. We 
conclude that the model ^ cannot explain the apparent 
discrepancy between energy scales in the optical spectrum 
and the conductivity measurements for the Bechgaard salt 
(TMTSF)2PF6. 

Parallel to the stacks of organic molecules the optical 
properties of (TMTTF)2X salts are clearly those of semi- 
conductors [T^HTl] . The optical spectrum displays sev- 
eral strong absorption features which are attributed to 
the coupling of electronic excitations with lattice vibra- 
tions. These transitions have a higher energy than the 
crossover energy above which excitations are effectively 
confined to a single stack and thus can be described by 
a one-dimensional model. A remarkable property of the 
(TMTTF)2X salts is that the optical gap is smaller than 
the Mott gap determined by photoemission experiments. 
For (TMTTF)2PF6 one observes that the strongest struc- 
ture in the optical spectrum is a relatively sharp peak at 
an energy (about 100 meV) significantly smaller than the 
Mott gap (about 200 meV). Obviously, this can be inter- 
preted as the signature of an excitonic transition below 
the gap for charged excitations. Thus, experimental ob- 
servations are compatible with the theoretical prediction 
of excitons in (TMTTF)2PF6 based on the model ^ and 
Mila's estimation for the appropriate parameters (es- 
pecially, the strength of the nearest-neighbor interaction 
V). Quantitatively, we find Ec ~ 160 meV and an exciton 
energy of toi w 60 meV using the parameters in Table d 
These energies are in very satisfactory agreement with the 
experimental values. Therefore, we conclude that excitons 
are present in the optical spectrum of (TMTTF)2PF6 [and 
probably other (TMTTF)2X salts] and explain the ob- 
servation of absorption features below the gap measured 
with photoemission spectroscopy. It would be interesting 
to have a direct experimental evidence for the presence 
of excitons in those salts. For instance, one could investi- 
gate the electro-absorption spectrum to demonstrate the 
presence of excitons as it was done for another quasi-one- 
dimensional material, polydiacetylene |32| . 

In conclusion, we have investigated the optical con- 
ductivity of the one-dimensional dimerized extended Hub- 
bard model at quarter filling using the dynamical density- 
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matrix rcnormalization group. Wc have found that the 
dimerization amphtude and the nearest-neighbor repul- 
sion (if strong enough to form excitons) determine the 
main features of the optical spectrum. Besides its influ- 
ence on the energy scale set by the Mott gap, the on- 
site repulsion plays a minor role only. Our study shows 
that this model cannot explain the optical spectrum in 
the Bechgaard salts (TMTSF)2X. It also shows that ex- 
citons probably contribute to the optical spectrum in the 
(TMTTF)2X sahs. 
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